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Constrained and stochastic variational principles for 
dissipative equations with advected quantities 

Xin Chen 1 , Ana Bela Cruzeiro 2 , Tudor S. Ratiu 3 


Abstract 

We develop symmetry reduction for material stochastic Lagrangian systems 
with advected quantities whose configuration space is a Lie group, obtaining de¬ 
terministic constrained variational principles and dissipative equations of motion 
in spatial representation. We discuss in detail the situation when the group in 
the general theory is a group of diffeomorphisms and derive, as an application, 
an MHD system for viscous compressible fluids. 


1 Introduction 

The goal of this paper is to develop a Lagrangian symmetry reduction process for a large 
class of stochastic systems with advected parameters. The general theory, which yields 
deterministic constrained variational principles and deterministic reduced equations of 
motion, is developed for finite dimensional systems. The resulting abstract equations 
then serve as a template for the study of infinite dimensional stochastic systems, for 
which the rigorous analysis has to be carried out separately. The example of dissipative 
compressible magnetohydrodynamics equations is treated in detail. 

The dynamics of many conservative physical systems can be described geometrically 
taking advantage of the intrinsic symmetries in their material description. These sym¬ 
metries induce Noether conserved quantities and allow for the elimination of unknowns, 
producing an equivalent system consisting of new equations of motion in spaces with 
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less variables and a non-autonomous ordinary differential equation, called the recon¬ 
struction equation. This geometric procedure is known as reduction, a method that 
is ubiquitous in symplectic, Poisson, and Dirac geometry and has wide applications 
in theoretical physics, quantum and continuum mechanics, control theory, and various 
branches of engineering. For example, in continuum mechanics, the passage from the 
material (Lagrangian) to the spatial (Eulerian) or convective (body) description is a 
reduction procedure. Of course, depending on the problem, one of the three represen¬ 
tations may be preferable. However, it is often the case, that insight from the other two 
representations, although apparently more intricate, leads to a deeper understanding 
of the physical phenomenon under consideration and is useful in the description of the 
dynamics. 

A simple example in which the three descriptions are useful and serve different 
purposes is free rigid body dynamics (e.g., [62, Section 15]). If one is interested in 
the motion of the attitude matrix, the material picture is appropriate. The classical 
free rigid body dynamics result, obtained by applying Hamilton’s standard variational 
principle on the tangent bundle of the proper rotation group 50(3), states that the 
attitude matrix describes a geodesic of a left invariant Riemannian metric on 50(3), 
characterized by the mass distribution of the body. However, as shown already by 
Euler, the equations of motion simplify considerably in the convective (or body) picture 
because the total energy of the rotating body, which in this case is just kinetic energy, 
is invariant under left translations. The convective description takes place on the Lie 
algebra so(3) of 50(3) and is given by the classical Euler equations for a free rigid 
body. Finally, the spatial description comes into play, because the spatial angular 
momentum is conserved during the motion and is hence used in the description of the 
rigid body motion. 

The present paper uses exclusively Lagrangian mechanics, where variational princi¬ 
ples play a fundamental role since they produce the equations of motion. In continuum 
mechanics, the variational principle used in the material description is the standard 
Hamilton principle producing curves in the configuration space of the problem that 
are critical points of the action functional. However, in the spatial and convective 
representations, if the configuration space of the problem is a Lie group, the induced 
variational principle requires the use of constrained variations, a fundamental result of 
Poincare [69]; the resulting equations of motion are called today the Euler-Poincare 
equations ([63], [62l Section 13.5], [H]). These equations have been vastly extended to 
include the motion of advected quantities ([131HS ;I6]) as well as affine ([36]) and non- 
commutative versions thereof that naturally appear in models of complex materials 
with internal structure ( [371 [39] ) and whose geometric description has led to the solu¬ 
tion of a long-standing controversy in the nematodynamics of liquid crystals ([Til Pf2]). 
The Euler-Poincare equations have also very important generalizations to problems 
whose configuration space is an arbitrary manifold and the Lagrangian is invariant 
under a Lie group action ([15]) as well as its extension to higher order Lagrangians 
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( [321 [28 i [29]). The Lagrange-Poincare equations turn out to model the motion of spin 
systems (|3T]), long molecules ( [2T], [30] ), free boundary fluids and elastic bodies ([34]), 
as well as charged and Yang-Mills fluids ([39]). There is also a Lagrange-Poincare the¬ 
ory for held theory ([38 : . :22]) and non-holonomic systems ([IS]). The Lagrange-Poincare 
equations also have interesting applications to Riemannian cubics and splines ([67]). 
the representation of images ([9]), certain classes of textures in condensed matter ([35] ). 
and some control ([30]) and optimization (|33j) problems. 

Variational principles play an important role in the design of structure preserving 
numerical algorithms. One discretizes both spatially and temporally such that the 
symmetry structure of the problem is preserved. Integrators based on a discrete version 
of Hamilton’s principle are called variational integrators ([64]). The resulting equations 
of motion are the discrete Eulcr-Lagrange equations and the associated algorithm for 
classical conservative systems is both symplectic as well as momentum-preserving and 
manifests very good long time energy behavior; see [olil loTj for additional information. 
There are versions of such variational integrators for certain forced ([IB]), controlled 
([68]). constrained holonomic ([581 Eo]), non-holonomic ([39]), non-smooth ([261 fl9j). 
multiscale m esi), and stochastic (0) systems. In the presence of symmetry, these 
systems can be reduced. However, today a general theory of discrete reduction in all 
of these cases is still missing and is currently being developed. If the configuration 
space is a Lie group, the first discretization of symmetric Lagrangian systems appears 
in [66], motivated by problems in complete integrability; for an in-depth analysis of 
such problems see [72] , 

All the above mentioned systems, both in the smooth and discrete versions, should 
have various stochastic analogues, depending on what phenomenon is modeled. The 
basic idea is to start with variational principles, motivated by Feynman’s path integral 
approach to quantum mechanics and also by stochastic optimal control. The latter 
has its origins in the foundational work of Bismut ([HE]) in the late seventies and 
in recent developments by Lazaro-Cann and Ortega ( [52l 531EHEB]). Non-holonomic 
systems have been studied in the same spirit in [33]. All of this work investigated 
mainly stochastic perturbations of Hamiltonian systems. A very recent approach on 
the Lagrangian side, in Euler-Poincare form, has been developed in [33], where both the 
position and the momentum of the system are (independently) randomly perturbed, 
as well as the Lagrangian. 

The present paper takes a different approach and is situated in the perspective 
of Feynman’s approach to quantum mechanics. Our Lagrangian is essentially (in the 
present paper, exactly) the classical one, but computed along stochastic paths, and 
their mean derivatives or drifts, in the spirit of m, he m ■ We derive deterministic 
equations for the drifts of these Lagrangian processes which satisfy now (dissipative 
and deterministic) partial differential equations whose second order terms correspond 
to the generator of the martingale part of the underlying paths. 

To illustrate our approach, let us consider a very simple example in the Euclidean 
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space M n . If g(t) denote deterministic integral curves with velocity field u, we have 
dg(t,x) = u(t, g(t,x))dt, g(0,x) = x. When the paths g are perturbed by a Brownian 
motion aW u (t) (a a constant), then they are no longer differentiable but we can still 
write dg u (t,x) = adW(t) + u(t, g u (t, x))dt, g u (Q,x) = x, valid almost everywhere in ui, 
as long as we interpret the differential d in sense of the Ito calculus. Heuristically, the 
increments g u (t + e,x) — g u (t,x ) have normal distribution with mean eu(t, g^t, x)) and 
variance ea 2 . Dynamically, the paths g follow “in the mean” the directions determined 
by the vector field u (the drift) but are subjected to a dispersion due to the Brownian 
motion effect. More precisely, we have, almost-everywhere in u, 

u(t, g u (t, x)) = lim - (E t g u (t + e, x) - g u (t, x)) 

£->■0 e 

where E t denotes conditional expectation given the past at time t. So, u still corre¬ 
sponds to a (generalized) derivative. The generator of the process g u (t) is the operator 
defined by Lf = ^-A/ + {u ■ V)/, i.e., we have 

E[/(0w(f,aO)] = f( x ) +E / [Lf(g u (s,x))]ds, 

Jo 

for every (compactly supported) smooth function /, where E denotes the expectation 
operator on the probability space whose points are aThe Laplacian term accounts 
for the dispersion. In this sense, L deforms the time derivative along the classical flow 
with velocity field u. 

Using such stochastic processes as “Lagrangian paths” one can derive variational 
principles extending those of classical mechanics. The critical points of action func¬ 
tionals, involving the above mentioned Lagrangians, solve almost surely the equations 
of motion. Those reduce to classical (Euler-Lagrange) equations in the limit a = 0, by 
construction. It follows, in particular, that as in this classical limit case, the Lagrangian 
encodes most of the qualitative properties of the (stochastic) system. 

In infinite dimensions, the Navier-Stokes equation was obtained in this manner, in 
the usual Eulerian (spatial) description ([2j l3l ITT]), via a stochastic Euler-Poincare re¬ 
duction on the group of volume preserving diffeomorphisms. The Lagrangian variables 
correspond to semimartingales. In this paper, we derive Euler-Poincare equations for 
stochastic processes defined on semidirect product Lie algebras and give the associ¬ 
ated deterministic constrained variational principle. In other words, we develop the 
semidirect Euler-Poincare reduction for a large class of stochastic systems. 

Plan of the paper. In Section [21 we recall some basic probability notions neces¬ 
sary for the rest of the paper and give the definition of the generalized derivative for 
topological group valued semimartingales. Section [3] contains the first main result of 
the paper, namely the stochastic semidirect product Euler-Poincare reduction for finite 
dimensional Lie groups, both in left and right-invariant versions. We give the deter¬ 
ministic variational principle and the reduced equations of motion. Section H] presents 
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the second main result of the paper, the reduction from the material to the spatial 
representation of the stochastic compressible magnetohydrodynamics equations. The 
stochastic reduction process recovers the standard deterministic equations in Eulerian 
representation. 

2 The derivative for semimartingales 

In this section, as in [3], we introduce the notion of generalized derivative for semi¬ 
martingales taking values on some topological groups. As we shall see later, the same 
definition holds if the semimartingale has values in some infinite dimensional group, 
such as the diffeomorphism group. 

2.1 Some probability notions 

We review in this subsection some basic notions of stochastic analysis on Euclidean 
spaces. We recall the concepts omitting the proofs, which can be found, for example 
in [47]. 

We denote R + := [0, oo[. Let (0, V, P ) be a probability space. Suppose we are given 
a family ( Vt)t&+ of sub-cr-algebras of V which is non-decreasing (namely, V s C Vt for 
0 < s < t) and right-continuous, i.e., D e >oVt+ e = Vt for all t G R + . We then say that 
the probability space is endowed with a non-decreasing filtration (Vt)tm+- 

A stochastic process X : R + x 12 —> R is (Vt)-adapted if X(t) is Vt- measurable for 
every t. Typically, hltrations describe the past history of a process: one starts with 
a process X and defines Vt n e >ocr{X(s), 0 < s < t}. Then the process X will be 
automatically (77)-adapted. 

A stochastic process M : R + x — y R is a martingale if 

(i) E\M(t,uj)\ < oo for all t > 0; 

(ii) M(t,u ) is (P f )-adapted; 

(iii) E s (M(t,u))) = M(s,uj) a.s. for all 0 < s < t. 

In the above definition, E denotes the expectation of the random variable, whereas 
E s (M(t,u>)) := K(M(t,u>)\V s ), for each s > 0, is the conditional expectation of the ran¬ 
dom variable M u (t), t > s, relative to the a—algebra (V s ), i.e., Id 3 u> i —> E s (M UJ (t)) e R 
is a P.,-measurable function satisfying 



Thus, condition (iii) is equivalent to E ((M u (t) — M u (s))xa) = 0 for all A G V s and all 
1,sGM satisfying t > s > 0, where xa is the characteristic function of the set A. 
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From now on we will often drop the probability variable co in the notations. In 
this paper we shall only consider processes defined in a compact time intervals [0, T] 
which are continuous (i.e., continuous with respect to the time variable t for almost all 
co G If). 

If a martingale M is continuous and E(M UJ (t) 2 ) < oo for all t > 0, we say that 
M has a quadratic variation {[M, Mj t , t G [0, T]} if M 2 (t ) — [M, M] t is a martingale, 
and [M, M] t is a continuous, adapted, a.s. non-decreasing process with JM, M] 0 = 
0. Such a process is unique and coincides with the following limit (convergence in 
probability), 

lim V] (M(t i+ i) — M(ti)) 2 

n—>-oo z —* 

where a n is a partition of the interval [0,t] and the mesh converges to zero as n —)■ oo. 

Actually, the definition of the quadratic variation requires only right-continuity 
of M. Moreover, for two martingales M and N, under the same assumptions and 
conventions as given above, one can also define their covariation 

lM,Nj t := lim V (M(t i+1 ) - M(t i ))(N(t i+1 ) - N(U)), 

n—>• oo zJ 

ti 5 ^ 1+1 

which extends the notion of quadratic variation. Clearly, 

2[M, Nj t = [M + N,M + Nj t - [M, Mj t - [IV, Nj t . 

More generally, one can consider local martingales. A stopping time is a random 
variable r : —y M + such that for all t > 0 the set {a; G | t(oo) < f} G V t . 
Then a stochastic process M is a local martingale if there exists a sequence of stopping 
times {r n ,n > 1}, such that liin^oo r n = oo a.s., and M n (t ) := M(t A r n ) is a 
square integrable martingale for all n > 1, where t A r n := min(f, r n ). We then define 
[M,M] t := [M",M"] t iff <r n . 

A real-valued Brownian motion is a continuous martingale W(t), t G [0, T], such 
that W 2 [t) — t is a martingale; or, equivalently, such that [W, W\t = t. 

A stochastic process X : Q x [0, T] —^ M is a semimartingale if, for every t > 0, it 
can be decomposed into a sum 

X{t)=X(0) + M{t)+A{t), 

where M is a local martingale with M(0) = 0 and A is a cadlag adapted process of 
locally bounded variation with A(0) = 0 a.s. (the definition requires that A be right- 
continuous with left limits at each t > 0; however, we consider only processes that are 
continuous in time, which is a standard assumption throughout this paper). 

More generally, M can be only a local martingale and A a locally bounded variation 
process. The definition of a local semimartingale follows the one given above for local 
martingales. 
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For a (local) semimartingale we define [X, X] t := [M, M\ t - 

Martingales and, in particular, Brownian motion, are not (a.s.) differentiable in 
time (unless they are constant); therefore one cannot integrate with respect to martin¬ 
gales as one does with respect to functions of bounded variation. We recall the defini¬ 
tion of the two most commonly used stochastic integrals, the Ito and the Stratonovich 
integrals. 

If X and Y are continuous real-valued semimartingales such that 
eQT \X(t)\ 2 dt +\Y(t)\ 2 d?j < oo, 


the ltd stochastic integral in the time interval [0, i], 0 < t < T, with respect to Y is 
defined as the limit in probability (if the limit exists) of the sums 


X(s)dY(s) 


lim V XiUKYiti+J-Yfa)) 

n—>oo L ^ 
ti j^i+1 


where a n is a partition of the interval [0,f] with mesh converging to zero as n —> 
oo. Moreover, when Y is a martingale such that E ^J Q T |A"(t)| 2 d[Y, < oo , then 

X(s)dY(s), t e [0, T], is also a martingale. 

The Stratonovich stochastic integral is defined by 

fx(s)SY(s)= lim £ (X(t i )+X(t H .i)) (y( 

/ n n—too i —* Z 

u k,t i+ i£cr„ 

whenever such limit exists. 

These integrals do not coincide in general, even though X is a continuous process, 
due to the lack of differentiability of the paths of Y. The Ito and the Stratonovich 
integrals are related by 

(2.1) f X S 5Y S = f X s dY s + \ f d{X, Yj s 

Jo Jo ^ Jo 

If / G C' 2 (1R), Ito’s formula states that 

(2.2) f(X(t)) = f(X(0)) + f f'(X(s))dX(s) + \f f"(X(s))dlX, X}, 

Jo * Jo 

This formula, for Stratonovich integrals, reads, 

f(X(t)) = f(X( 0))+ f f'(X(s))5X(s) 

Jo 
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One advantage of Stratonovich integrals is that they allow the use of the same rules 
as those of the standard deterministic differential calculus. On the other hand an 
Ito integral with respect to a martingale M is again a martingale, a very important 
property. For example, we have, as an immediate consequence, that E s f X(r)dM (r) = 
0 for all 0 < s < t. 

In higher dimensions, the difference between the Stratonovich and the Ito integral 
in Ito’s formula is given in terms of the Hessian of / (see subsection 12.2p . In fact, 
suppose that X be a Revalued semimartingale; then Ito’s formula in d-dimensions 
(see also (12.2(1 ) states that, for every / G C 2 (M d ), 

f(x(t)) = f(x(0)) + V f t aj(X{s))dx i {s) + lY' f 

(2.3) = f(X( 0)) + J2 f dif(X( S ))t>X\ s ) 

i =1 

For Brownian motions W\ i = 1... ,k, we have 

(2.4) d\W\ W% = Sijdt 

where 8^ denotes the Kronecker delta symbol, and, as the covariation of semimartin¬ 
gales is determined by their martingale parts, the following identities hold (see, e.g., 

El). 

(2.5) d\W\t J, = 0 Vi = 1,.... d, = 0, 

where t(t) = t is the identity (deterministic) function. 

2.2 The generalized derivative for (topological) group valued 
semimart ingales 

Let G denote a topological group, endowed with a Banach manifold structure (possibly 
infinite dimensional) whose underlying topology is the given one, such that all left (or 
right) translations L g (resp. R g ) by arbitrary g 6 G are smooth maps; L g h := gh, 
R g h := hg, for all g,h G G. Given v G T e G, denote by v L (resp. v R ) the left 
(resp. right) invariant vector held whose value at the neutral element e of G is v, i.e., 
v L (g) \= T e L g v (resp. v R (g) := T e R g v ), where T e L c ■ T e G —> T g G is the tangent map 
(derivative) of L g (and similarly for R g ). The operation [ni,n 2 ] := (e), for any 

Ui,t >2 £ T e G, defines a (left) Lie bracket on T e G. When working with right invariant 
vector fields, we shall still use, formally, the left Lie bracket defined above, i.e., we 
shall never work with right Lie algebras; the bracket defined by right invariant vector 





fields is equal to the negative of the left Lie bracket defined above. Denote, as usual, 
by a,d u v := [u,v] the adjoint action of T e G on itself and by ad* :T*G —$■ T*G its dual 
map (the coadjoint action of T e G on its dual T*G ). 

Suppose that V is a left invariant linear connection on G, i.e., \7 v lv 2 is a left 
invariant vector field, for any vi,v 2 G T e G. Then we define V Vl v 2 : = V^uf^e) for all 
V \, v 2 G T e G. If right translation is smooth, in all the definitions above, we can replace 
left by right translation in a similar way. The left invariant connection V is torsion 
free if and only if 


V, ;i v 2 - S7 V2 v i = [ui, v 2 \ , for all v±,v 2 G T e G 

as an easy verification shows. For a fixed g± G G, let T 92 L gi : T 92 G — > T gi92 G be the 
tangent map (or derivative) of L gi at the point g 2 G G. 

Let G be endowed with a left invariant linear torsion free connection V. The 
corresponding Hessian Hess f(g) : T g G x T g G — > M of / G C 2 (G ) at g G G is defined 
by 

(2.6) Hessf(g) (ui, n 2 ) := Viv 2 f(g) - v 1: v 2 G T g G, 

where h,, i — 1,2, are arbitrary smooth vector fields on G such that Vi(g) = n,. Since 
the connection is torsion free, Hess f(g) is a symmetric M-bilinear form on each T g G. 
In addition, Hess/ = V 2 / = Vd/ (see, e.g., [Tj or [23]) is the covariant derivative 
associated with V of the one-form d/, where d denotes the exterior differential. 

A semimartingale with values in G is a stochastic process g : Q x M + —> G such 
that, for every function / G C 2 (G), fog : H x M + —)• R is a real-valued semimartingale, 
as introduced in subsection 12.11 (see, e.g., |23| for the case of finite dimensional Lie 
groups). 

A semimartingale with values in G is a V- (local) martingale if 

t 1 —» - f(g w ( 0)) - Hess f(g u (s))dlg u , g^] ds 

is a real-valued (local) martingale, where | gu,g u \t is the quadratic variation of g u . If 
G is a finite dimensional Lie group, then we have the following expression 

/ p; 1 ^), / p; 1 ^)! , 

Jo Jo JJ t 

where P t : T gu] ^G —$■ T g ^G is the (stochastic) parallel translation along the (stochas¬ 
tic) curve 1 1 —y gu{t) associated with the connection V, see e.g., [23] or [4T] . Moreover, 
for some infinite dimensional groups G (for example the diffeomorphism group on a 


d\g u ,g u \t := d 
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torus), the quadratic variation is also well defined: we refer the readers to [3] jT7] for 
details (see also Section 4 of this paper). 

For a G -valued semimartingale g u (-), suppose there is a, possibly random, process 
: (t,u) H * TG, such that v(i) G T gut ^G a.s., and for every / G C°°(G), the 

process 

N[ := f(g u {t)) -/(&,(()))) - Hess/(^(s)))d[^,^] s - (d f(g u (s)),v u (s)) ds 


is a real-valued local martingale. Then we define the V -generalized derivative of g u {t) 
by 


(2.7) 


dt 


■= V w (t). 


We remark that by (12.71) . the V-generalized derivative is well defined for semimartin¬ 
gales with values in a finite dimensional Lie group as well as in some infinite dimensional 
groups (the diffeomorphism group on a torus for example), see e.g. |3j. Moreover, when 
G is a finite dimensional compact Lie group, we have the following equality (see [23] ) 


( 2 . 8 ) 


dt 


■ =Pt 


lim E t 

e^O 


Vw(t + e) 
e 



e T gu,(t)G , 


where 

v u {t) = I P T'SgM e T e G. 

Jo 

The conditional expectation E t plays a major role in the definition of the V-generalized 
derivative, since it eliminates the martingale part of the semimartingale. Therefore, 
the velocities are given by the drift (the bounded variation part) and the diffusion 
part (the martingale) can be seen a stochastic perturbation; in other words, the drift 
determines the directions where the particles flow, the martingale part describes their 
random fluctuations. 

Moreover, according to the definition, if a G -valued semimartingale g u (t) satisfies 
9 9 d ^ = 0, then g u (t) is a V-martingale. 

In Euclidean space, endowed with the standard Levi-Civita connection of the con¬ 
stant Riemannian metric given by the inner product, (12.7[) coincides with the definition 
of the generalized derivative for the Euclidean valued semimartingale considered in El 
and [74] , [75]. 

Given Revalued Brownian motion W w {t) = (W^(t), ..., W%(t)), t G [0,T], vectors 
Hi G T e G, 1 < i < k, and curve u G G 1 ([0, T]; T e G), we consider the following 
Stratonovich SDE on G, 


dg u {t) = T e L gu{t) (eLi - | V Hi Hidt + u(t)dtj , 

9u{0) = e; 
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here, for each v,w G T e G , we recall that V„w G T e G by X7 v w := (Vviv) (e), where w 
and w are the left invariant vector fields (or right invariant, if the connection is right 
invariant) on G, such that v(e) = v, w(e) = w. 

Since dl(T e L gu ( t )Hi),W„(t)}t = T e L gui (t)(V Hi Hi)dt, formula (T2TT1) shows that equa¬ 
tion (12.911 is equivalent to 

(2 10) f dgM = T e L„ At) (El, HidW&t) + u(t)dt ) , 

\ &,(0) = e. 

If G is a finite dimensional Lie group, there exists a unique strong solution for (12.9ft 
(c.f. 07|, H). When G is the diffeomorphism group on a torus a weak solution still 
exists ([3], [17]). 

Applying Ito’s formula for G-valued semimartingales (see [23] for the case where 
G is finite dimensional and |3J, section 4.2 for the case where G is the diffeomorphism 
group on a torus), for every / G C 2 (G) we have, 

/(&>(*)) = /(&;(0)) + N( + ]- [ Hess f(g u (s))dlg u ,g u } s + I T e L gUs) (u(s))f(g u {s))ds 

* Jo Jo 

where N{ is a martingale. Actually, this last equality, valid for each / G C 2 (G), is a 
characterization of the solution of the stochastic differential equation (12.9[) (or (12.101) 1. 
in a weak sense. 

Clearly, when G is a finite dimensional Lie group, by definition (12. 7|) . we have 
(2-11) ^ — T e L gu{t) u{t). 

3 Stochastic semidirect product Euler-Poincare re¬ 
duction 

In this section, inspired by [3], we generalize the deterministic semidirect product Euler- 
Poincare reduction, formulated and developed in in [45] . to the stochastic setting. 

Left invariant version. Let U be a vector space, U* its dual, and ('> ')u : U*xU -G R 
the duality pairing. Let G be a group endowed with a manifold structure making it 
into a topological group whose left translation is smooth. As discussed in subsection 
12.21 the tangent space T e G to G at the identity element e G G is a Lie algebra. 
Assume that G has a left representation on U; therefore, there are naturally induced 
left representations of the group G and the Lie algebra T e G on U and U*. All actions 
will be denoted by concatenation. Let (•, -) T G : T*G x T e G —* M be the duality pairing 
between T*G and T e G. Define the operator o : U x U* * T*G by 

(3.1) [a o ct, v) t ^q '.= — {va, cl) u = (ct, va) v , v G T e G, a G U, aEU*. 
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In fact, a o a is the value at (a, a) of the momentum map U x U* —>■ T*G of the 
cotangent lifted action induced by the left representation of G on U. 

Let 5?(G) denote the collection of G- valued continuous semimartingales defined on 
the time interval t G [0, T]. Given a (left invariant) linear connection V on G, a point 
«o G U*, and a (Lagrangian) function l : T e G x U* —> M, dehne the action functional 
JV’“o,i . y(G) x y{G) M+ by 

(3.2) J v '“«’'(9i(-),£(-)) := E £l (T s ,j t) L gUtr ,^^-,am'j dt, 

where e ^(G) and 

(3.3) a(t) :=E[a w (t)] G U*, cfo(f) :=sS(t) _1 ao- 

For every (deterministic) curve ^(-) G C' 1 ([0,1]; T e G) satisfying ^(0) = ^(T) = 0 
and £ G [0,1), let e £ ^(-) G G 1 ([0,T];G) be the unique solution of the (deterministic) 
time-dependent differential equation on G 


(3.4) 


— £ T e L es ^(t)£(t) , 
e £ii? (0) = e. 


Note that this system implies eo^ (i) = e for all t G [0, T], 

We say that (fli(-), fl£(')) e -^(G) x y{G) is a critical point of J V,Q W if for every 
(deterministic) curve g{-) G C ,1 ([0, l];T e G) with ^(0) = g(T) = 0, we have 


(3.5) 

where 

(3.6) 


d 

de 



= 0 


9l,e,M := 9l(t)e e>8 (t), t G [0, T], 


1,2, £G [0,1). 


We emphasize the particular form of these deformations in the Lie group: they corre¬ 
spond to developments along deterministic directions p{t). 


Remark 3.1. As will be seen in the applications in the next section, the reason we 
choose two different semimartingales (gd(-), dt{')) i n the variational principle (13. 5 ft is 
that the viscosity constants in the equations for the variables in T e G and U* may be 
different. <0 


From now on in this section we assume that both G and U are finite dimensional. 
Theorem 3.2 still holds for some infinite dimensional models, see section 4 for the 
case of the diffeomorphism group on the torus. We fix HJ 1 G T e G, with m = 1,2, 
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j = 1 as well as W™(t) = W™’ km (t)), m = 1,2 independent 

valued Brownian motions (notice that, in the infinite dimensional case, we may not be 
able to fix a priori the Brownian motions in order to solve the corresponding stochastic 
differential equations). 

We consider semimartingales g™(-) G y{G) 1 m — 1,2, defined by 


i£(t) = T e L, s(f) (E-:, + u(t)dt) , 

C(0) = 


for some w(-) G C 1 ([0,T];T e G). Note that u(-) is not given (and is the same for m = 
1,2); we shall see below that it is the solution of a certain equation when (c/^(-), <?£(•)) 
is critical for J v,Qo,z . 

In the theorem below we use the functional derivative notation. Let V be (a possibly 
infinite dimensional) vector space and V* a space in weak duality (•, •) : V* x V —> R 
with V ; in finite dimensions, V* is the usual dual vector space, but in infinite dimensions 
it rarely is the topological dual. If / : V —> M is a smooth function, then the functional 
derivative G V*, if it exists, is defined by lim f ^ 0 Tu+d)-.fuo = (|{, 6) for all a,b E V. 
Thus, in the theorem below G TfG and G U are the two partial functional 
derivatives of / : T e G x U* —> M. 

Theorem 3.2. Assume G is a finite dimensional Lie group endowed with a left invari¬ 
ant linear connection V. Let U be a finite dimensional left G-representation space and 
l : T e G xl/WKfl smooth function. Suppose that the semimartingales g™(-) G 5^{G), 
m = 1, 2, have the form (13.711 and their deformations are of the form (13.61) . 

(i) Then (<?i(-), 9u(')) a critical point of J v ’ a °’ 1 ( see (13.2jl ) if and only if the ( non- 
random ) curve u(-) G C 1 ([0, T]; T e G) coupled with «(•) G C' 1 ([0, T]; U*) (defined 
by (13.311 ) satisfies the following semidirect product Euler-Poincare equation for 
stochastic particle paths: 

('3 8') / = ad Si£ + £ oa + K (£) > 

1 = I Ejii H j (#/«(*)) ~ u 2 (t)a(t), 

where G T*G and G U denote the two partial functional derivatives of l, 

the operation o is defined by (13.ljl . a(-) by (13.3ft . 



j= 1 
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(3.9) 








and the operator K : T*G —* T*G is defined for all /i G TfG and v G T e G by 


(3.10) 


(A'(fi).ti) 



(ii) The first equation in (13.81) is equivalent to the dissipative Euler-Poincare varia¬ 
tional principle 

(3.11) 5 f l(u(t), a(t))dt = 0 

Jo 

on T e G x U* defined by l : T e G x U* —> M, for variations of the form 


(3.12) 


{ 5u = v + [u\ u] + K*(v ) = v + [u, v\ 
5a = —va 


1 

2 


ki 




■ 3 = 1 


+ K*(v), 


where v[t ) is an arbitrary curve in T e G vanishing at the endpoints, i.e., v (0) = 0, 
v{T) = 0. Note that this variational principle is constrained and deterministic. 

Proof, (i) Step 1. We start by proving prove that a{t) satisfies the second equation 
in (13.81) . We have 


d (&(*)) 1 = - T eR(gi(t))-^T g 2 it) L {g 2 {t)) -idgl{t) 
and, replacing dg^{t) by its expression, 


(3.13) 


«(<))-' = T e R ( ,lm-' (E‘i, (-HfSWp(t) + \V„iH]di) - u(t)dt ) , 

m = e > 


We now derive the stochastic differential equation satisfied by a u (t) := gf{t) 1 a 0 : 


(3.14) da w (t) = d (gl(t) ^o) = [-T g 2 {t) L g 2 {t) -idgl(t)\ g£(t) x a 0 

= - jt (h? (gl(t)-'a „) swp(t) - (sHtr'a o) dtj 

- u (t) dt 


or, equivalently, 

(3.15) da w (t) = d (< ql(t)~ l a 0 ) 
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D rr, D 


fc 2 


X] ( - (//|a 0 ) SW^(t) + - a 0 dt ) - u(t)a 0 dt 

.1=i 


Since we assume U* to be a finite dimensional vector space, the difference between the 
Stratonovich and Ito integrals (see (12.11) 1 yields 


fc 2 


E ( H J (sK.tr'co)) w?(t) 


3 = 1 


k2 


= E (( H i „)) (t) + -d[y 2 (aJ(-)-‘ao). w^] t 

By the same procedure as in (13.141) . the (local) martingale part of H 2 (g^-)- 1 a 0 ) is 
equal to — Xu=i fo H 2 H 2 (g 2 (s)~ 1 a 0 )dW 2 ’ l (s). Therefore, by (12.41) and (12.5[) we derive 

E <il Hj (£(■)-'<* o) . Wy), = - E H](H]( g l{t)- l a a ))dt 

1=1 1=1 


Using (13.14j) we have, 


fc 2 

(3.16) dc^f) = - X] (iZ?(a w (i))dW^W) 
l=i 

1 fc 2 1 fc 2 

+ ^ H i ( H i (“«(*))) dt + 2 V H ^ R2j ) dt ~ u (*) (“«(*)) 

Z l=i Z l=i 


or, equivalently, 

(3.17) da w (t) = T e R g 2 p)-i 

+ T e R g i{t) 




E - (y 2 <*o) diry(i) 


-l=i 
fc 2 


-l=i 


y: ( —Hj(HjOio)dt + - (VffiHfj a 0 dtj - u(t)a 0 dt 


Note that (13.141) . respectively (I3.16p . are equations for a u (t) = (g^(t)) 1 a 0 , a new 
variable that is introduced by the reduction process. On the other hand, (13. 15j) and 
its Ito version (13. 17jl are equations for the right logarithmic derivative of g^(t) with 
known right hand side (once u(t) is determined). Thus, if one wants to solve for g^{t), 
this second group of equations are more convenient, even though they do not use the 
reduced variables. 
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Taking the expectation in (3.10), as Ito integrals have vanishing expectation values, 
we have 

(3.18) - E [fi 2 (i)(S„(«))] , 

3 = 1 

where u 2 (t) is defined in 03.9[) . Since the Lie algebra representation on U* commutes 
with integration over hi (as linear maps), we deduce that a(t) := E [fl£(£) _1 a!o] G U* 
satisfies the equation 


(3.19) 





— u 2 (t)a(t). 


Step 2. We prove the first equation in f!3.81) . Recall from (13.4p that, for every 
(deterministic) curve g(-) G C Yl ([0,1]; T e G) satisfying ^(0) = g(T) = 0 and e G [0,1), 
e £rf (-) G C 1 ([0, T] \ G) uniquely solves the following differential equation on G 

e e,^(0) — &■ 

Recall also the notation g l Uj£:£ (t) := g l u (t)e £ j(t), t G [0,T], i = 1,2. 

By (31 Lemma 3.1], we have 


(3.20) 


d_ 

de 


e erf (f) =0{t), 

£=0 


d_ 

de 


^8 


(t) 


-1 


£=0 


-8 if) 


Since this computation is important in the proof, for the sake of completeness, we 
recall it below. Denoting by jjj and jf the covariant derivatives, induced by V on 
G, along curves parametrized by t and e, respectively, since the torsion vanishes, the 
Gauss Lemma yields 


(3.21) 


D d D d D, .. x 

Dide ee * (t) = D~eJt e£ ^ t] = 

= T e L e e ^{t)8f) + £—(T e L €e s(f) g(t)) 


Taking e = 0, since e 0 ^(t) = e for all t, we obtain _ =0 e £ji9 (t) = g(t). How¬ 
ever, t H> ^|__ 0 e £i3 {t) is a curve in the vector space T e G and hence ^^| g _ 0 e £ ^(t) = 
Ml=o e *At) = 8f)- The equality in (13.201) is then a consequence of ^(0) = 0 
and o e £rf (0) = 0. Finally, since 


d 

de 


e e&{t) 


-1 


—Tp_R-i t .\T F „L- 


d 


eG(0 £ ’^ e G(0 de ' 8 


(t), 
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the second equality in (I3.20[) follows from the first. 

Next, we need the following identity from the proof of [3| Theorem 3.2] 


(3.22) 


d_ 

de 


e=0 


T 'i 

yu,e,g 


d)-^9i,e, s dy 


dt 


dgft) 

dt 


+ adui(t)^(t) 


n,l 

2 ( v a+ Vffi (ad 0 {t)H])^j , 

3 =1 


where u l {t ) is given by (13.91) . To prove it, introduce H)' £ (t) := Ad e -i^ Hj , use (13.7|) . 
Ito’s formula, and recall that e £ ^(t) is of bounded variation (the time derivative of a 
product obeys the usual Leibniz formula only if at least one of the factors is of bounded 
variation). We get 


fci 


1 


fci 


= £ T e L A , im H^(t)SW^(t) + T.L^ m Ad e - i(() u(i)) dt 

j =i V j=i 

+ r d-'eLe e , s (t)- 1 &E,g(d)dt 

k\ 

£ Ad ed(1) u\t)dt 

3 =1 

or, written in the Ito form, 


fci 




0=1 


/Cl 


+ ^2 ^ Hj’ £ (t)Hj ft) + Ad e -i( t ) u (£) + T e L ees (j^-ie £tB {t) j dtj . 

By Ito’s formula (for hnite dimensional manifolds) and by (13.4p . for every function 
/ G C 2 (G), the process Nf defined as follows is a real-valued local martingale, 


N ft,ft) -=f{9l,E, 3 ft)) - f{al,E, 3 m 


fc'i r t 


3=1 J0 


n£ / )(gl„Js))ds 


T « L 9l,„(*)( Ad ««-J(.) fi‘W + £>(»))/) ((ijs))ds 
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= f(9l,e, £ (t)) - f{9l, es { 0)) 


^ i=l 7 ° 


1 


fci r t 


E / T ^i ^ H] ^ s) H^(s)f(gl^s))ds 


j =i Jo 


T e L <e>)( Ad erj( S ) “ 1 ( S ) + £ i( S ))/) (^,^( S )) rfS 


Thus, formula (12.lip for the V-generalized derivative of the semimartiugale (cu,t) i—>• 
9 u{t)eg^) I s applicable and we obtain 




1 ft) 1 fc l 

E V„.,. (t) Hp(i) + Ad e -i (t) «>(«) + £*(S). 


dt 


2 — -n 
i=i 


To prove (I3.22j) . we take ^| £ _ 0 of both sides. The third term yields ^(t). To compute 
this derivative of the second term, use (I3.20p and e £ii? (0) = e to conclude 


d 

de 


Ad eri(t) “ 1 ( t ) = — ac W) u 1 (t) = ad^i {t)3 (t) 


e=0 


Finally, to compute this derivative of the first term, note that F/y’ u (f) = Hj for all t 
and that, again by (I3.20p . Hj' s (t) = — ad £ (t) Hj, which yields 


d 

de 


^- V ^ 1 1 ( ad ^(b H j 


£=0 


- 


thus proving (13.221) . 

Since gl^t) := gl(t)e eig (t) and e 0t3 {t) = e for all t e [0, T], we conclude gl^it) = 
gl,(t), for all t e [0, T], i — 1,2. Therefore, 


. , d 

(3 ’ 23 > s 


£=0 


^0 = -gl(t) 1 ^ 


£ = 0 


at^it) giit) a 0 


d 

de 


£=0 


e Crf (i)^ 1q o ™ -g(t)gl(t) W 


With (I3.22p and (I3.23j) . we can compute the e-derivative of J v ’ ao,/ 9u,e,#(') 

at e = 0. We get from (13.2[) 


d 

de 


J 


£=0 


^ 7,q °^ ffid ( •) o 2 

\Hu},E,£\ ) ; Voj,e 
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- J Q (j(t) ( E [slit) V)]) 

= ^ ^ + + ° + A " (^W) > ^(*)^ di > 


where the last equality is obtained integrating by parts the first summand (the bound¬ 
ary terms vanish since ^(0) = g(T) = 0), and applying the definitions of ad*, o (see 
f|3.ip ). and K (see (13.10P ). 

Since g E C' 1 ([0, T]; T e G) is arbitrary, (gd(-), <?£(■)) is a critical point of J V,Q °’ 1 if 
and only if u(-) E C' 1 ([0, T]; T e G ) satisfies the second equation in (13.81) . 

This proves statement (i). 

(ii)The expressions of the variations (13.121) have been already found in the previous 
computations. Using them in the variational principle (13.111) . we obtain 


5 / l(u(t),a(t))dt = 


lf3~T2l) 


lL M ) + ( Sa ’^)) dt 

—, u + ad s i v + K*(v) 
ou 


d (5l\ * 51 


va ’t) ]dt 

51 \ 51 

-\ + -« a ,v) d t. 


This vanishes for all smooth curves v(t) E T e G vanishing at the endpoints if and only 
if the first equation in (13.8p holds. □ 

Remark 3.3. As we shall see later, the conclusion of Theorem 13.21 still holds when G 
is the diffeomorphism group of a torus and the action of G on U* is the pull back map. 

If G is a finite dimensional Lie group and U a finite dimensional vector space, then 
a(-) and a(-) defined by (13.3p satisfy some SDE and ODE respectively. However, when 
G is the diffeomorphism group, as will be seen in the next section, a(-) and a(-) satisfy 
some SPDE and PDE, respectively. 0 

Instead of (13.7p . let us consider semimartingales h™(-) E S^{G), m — 1,2, defined 
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by 


dhZ(t) = T e L hsm (£‘ =1 , 

W) = e, 


for «(•) G C 1 ([0, T}\T e G). Accordingly, define 

(3.25) a(t) : = E[d^(f)] G U*, a u (t) := /i^,(f) _1 a 0 . 

Then, proceeding as in the proof of Theorem 13.21 we get the following result. 

Theorem 3.4. Under the conditions of Theorem 3.2, suppose that the semimartingales 
hff{-) G y(G), m — 1,2, have the form (13.24[) and their deformations are of the form 

(BSD- 

(i) Then (/&*(•), ^(-)) is a critical point of J v,ao ’ 1 ( see (13.21) 1 if and only if the ( non- 
random ) curve u(-) G C 1 ([0, T]; T e G) coupled with «(•) G C' 1 ([0, T\; U*) ( defined 
by (13.25R ) satisfies the following semidirect product Euler-Poincare equation for 
stochastic particle paths: 


(3.26) 


= l J2jii H j - u 2 (t)a(t), 


(ii) The first equation in (I3.20p is equivalent to the dissipative Euler-Poincare varia¬ 
tional principle 

(3.27) 5 f l(u(t), a(t))dt — 0 

Jo 

on T e G x U* defined by l : T e G x U* —> M, for variations of the form 

5u = v+ [u 1 , v] + K*(v) = v + [u, v] + K*(v), 
da = —va 


where v(t ) is an arbitrary curve in T e G vanishing at the endpoints, i.e., n(0) = 0, 
v(T) = 0. 


Right invariant version. Due to relative sign changes in the equations of motion 
and the dissipative constrained variational principle, with a view to applications for 
the spatial representation in continuum mechanics, we give below the right invariant 
version of Theorem 13.21 

Let G be a Lie group endowed with a right invariant covariant derivative V, i.e., 
V v rV 2 is a right invariant vector field, for any Vi,v 2 G T e G, where v[ { G X(G) is 
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the right invariant vector field defined by vf{g) := T e R g Vi , i = 1,2. Recall that 
[vi,V 2 ] R = — [vf, v^\, where on the left hand side the bracket is in g, hence defined 
via left invariant vector fields (the standard left Lie algebra of G), and on the right, it 
is the usual Jacobi-Lie bracket of smooth vector fields on G. 

Suppose that G acts on the right on a vector space U (and we will write the action 
of g G G on u G U by ug and similarly for the induced infinitesimal g-representation). 
In the procedure leading to Theorem 13.21 and its proof, we interchange all left actions 
and left translation operators by their right counterparts. Thus the action functional 

jV,a 0 ,l 

is defined in this right invariant case by 


(3.29) J v ’“»''( 9 y-).^(-))-=E / 


9 “ (t) glgle.no), 


where T g i^R g i^-i denotes the differential of the inverse of right translation R g i(p-i 
at the point g]j{t) G G and 


(3.30) 

a{t) \= E [a w (t)] G U*, atuit) := a 0 gl(t) 1 . 

As before, for every (deterministic) curve G C^QO, 1]; T e G) satisfying ^(0) = 

g(T) = 0 and £ G [0,1), let e £rf (-) G G 1 ([0,T];G) be the unique solution of the 
(deterministic) time-dependent differential equation on G 

(3.31) 

i dt e£ — £ TeRe e ,j,(t)jj{t), 

\ gJ 0 ') = e. 

Define 

:= e eig (t)gl(t); 


we can consider (right invariant) critical points of as before. Moreover, formulas 

(13.7[) (j3. 19j) still hold with all the left translations and left g-actions replaced by right 
ones. For example, the right version of the time evolution (13. 19(1 of a{t) (given this 
time by (13.30(1 )) takes the following form for right invariant linear connections on G and 
right G-representations 


(3.32) j t a(t) = E) - a(t)u 2 (t), 

3 = 1 

while the definition (13.9(1 of u 2 (t) remains unchanged. The semimartingale g 2 (-) G 
S^(G) is the solution of the following two equivalent stochastic differential equations 
(in Stratonovich and ltd form) 

(3.33) d&u(t) = d {a Q gl(t)~ l ) 
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k 2 / \ 

X] ( - (« 0 H?) 5W 2 J(t) + -a 0 dt - a 0 u(t)dt 

U=i ^ ' 

it - (“o h i) d < j ^) 

Li=i 


+ T e L gl(ty 


k 2 


J=1 


Y f 7 }(ctoHj)H 2 dt + ^a 0 (v dt ) - a 0 u(t)dt 


Finally, in the proof of the analogue of (13.22)1 (on the right hand side the second and 
third summand change sign), one uses H^' £ (t) := Ad eejj (t) Hj. 

The right invariant version of Theorem 13.21 is the following. 


Theorem 3.5. Assume G is a finite dimensional Lie group endowed with a right 
invariant linear connection V. Let U be a finite dimensional right G-representation 
space. Suppose that the semimartingales gff) £ m — 1,2, have the following 

form, 


(3.34) 


ACM = TeR g -(t) ( Eh, f) - \V„ T Hfdt) + u(t)dt ), 

9”(0) = e, 


and that their deformations are given by (13.6p . 

(i) Then (^i(-), 9u(')) a critical point of J v,ao ’ 1 given by (13.291) if and only if the ( non- 
random )I curve uf) e C 1 ([0, T]; T e G) coupled with «(•) G C 1 ([0, T]; U*) (defined by 
(13.30)1 ) satisfies the following semidirect product Euler-Poincare equation for stochastic 
particle paths: 


(3.35) 


|«M = I E‘i, ymj) H] - 


where cn(t), u 1 , and K are defined by fl3.30p . (13. 9p . and (13.101) . respectively. 

(ii) The first equation in (I3.35P is equivalent to the dissipative Euler-Poincare variational 
principle 


(3.36) 


5 f l(u(t),a(t))dt = 0 
Jo 


on T e G x U* 

defined by l : 

T e G xf/MM, for variations of the form 




- ki 

(3.37) 

5u = v — [ft 1 

,v\ - K*(v) = v + [u, v] - - 

J= 1 


^ 5a = —av 




K*(v), 
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where v(t) is an arbitrary curve in T e G vanishing at the endpoints, i.e., u(0) = 0, 
v(T) = 0. 

Remark 3.6. Under some special conditions, the operator K(u ) takes a familiar form. 
More precisely, if V is the (right invariant) Levi-Civita connection and if we assume 
that V H t Hi = 0 for each i, we have, 

(3.38) K(u) = ^ (V/TiVtfiU + R (u, Hf)Hf), \/u G 0 , 


where R is the Ricmannian curvature tensor with respect to V. In particular, if {//,} 
is an orthonormal basis of 0 , then K(u) = — |(A u + Ric(w)), where A u := AI/(x)| ff=e 
for the right invariant vector fields U(g) := T e R g u, \/u G g, g G G. 

A proof can be found in [3], Proposition 3.5. 0 

The right invariant version of Theorem 3.5 reads: 


Theorem 3.7. Assume G is a finite dimensional Lie group endowed with a right 
invariant linear connection V. Let U be a finite dimensional right G-representation 
space. Suppose that the semimartingales hff) 6 5?{G), m = 1 , 2 , have the following 
form, 


(3.39) 


dhZ(t) = T e R hw f £‘" (HfSW^(t))+u(t)dt), 
'C(O) = e, 


and that their deformations are given by (13.61) . 

(i) Then (hf(-), is a critical point of given by (j3.29j) if and only if the ( non- 

random ]I curve u(-) G U 1 ([0, T]; T e G) coupled with a(-) G U 1 ([0, T]; U*) (defined by 
(13.301) ) satisfies the following semidirect product Euler-Poincare equation for stochastic 
particle paths: 


(3.40) 





\a(t) = i Y .%i (“( 4 )y 2 ) H i - a(t)u 2 (t) 
where K is defined by fl3.10p and aft ) is defined by 


(3.41) a(t) := E [a^t)} G U*, a u (t) := a 0 h 2 u {t) 1 . 


(ii) The first equation in (13.40j) is equivalent to the dissipative Euler-Poincare variational 
principle 


(3.42) 


S / l(u(t), a[t))dt = 0 
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on T e G x U* defined by l : T e G x £/*—>• R, for variations of the form 


(3.43) 


5u = v — [it 1 , -y] — K*{y) = v + [w, u] — K*(v), 
5a = 


where v(t) is an arbitrary curve in T e G vanishing at the endpoints, i.e., u(0) = 0, 
v(T) = 0. 


4 Application to the MHD equation 

We begin by recalling from [20] and [M] the necessary standard facts about the group of 
diffeomorphisms on a smooth compact boundaryless n -dimensional manifold M. Then, 
when we present the compressible MHD equations, we shall specialize to the periodic 
case, i.e., we shall take M = T 3 , the usual three dimensional flat torus. Of course, all 
of the discussion below can be easily extended to smooth manifolds without boundary. 

Let M be a smooth compact boundaryless n-dimensional manifold. Define 

G s := {g : M —> M is a bijection | g,g~ l E H S (M, M)} , 

where H S (M , M) denotes the manifold of Sobolev maps of class s > 1 + ^ from M to 
itself. The condition s > | suffices to ensure the manifold structure of H S (M, M); only 
for such regularity class does the notion of an H s - map from M to itself make intrinsic 
sense. If s > 1 + | (the additional regularity is needed in order to ensure that all 
elements of G s are C 1 and hence the inverse function theorem is applicable), then G s 
is an open subset in H S (M, M ), so it is a C°° Hilbert manifold. Moreover it is a group 
under composition between maps, right translation by any element is smooth, left 
translation and inversion are only continuous, and G s is a topological group (relative 
to the underlying manifold topology) (see [20], [61]); thus, G s is not a Lie group. Since 
G s is an open subset of H S (M , M), the tangent space T e G s to the identity e : M —> M 
coincides with the Hilbert space X S (M) of H s vector fields on M. Denote by X(M) the 
Lie algebra of C°° vector fields on M. The failure of G s to be a Lie group is mirrored 
by the fact that X s (M) is not a Lie algebra: the usual Jacobi-Lie bracket of vector 
fields, i.e., [X, Y}[f] = X[Y[f]}-Y[X[f}} for any X,Y E X(M) and / E C°°(M), where 
X[f] := d f(X) is the differential of / in the direction X , loses a derivative for finite 
differentiability class of vector fields and thus [■, ■] : X s (M) x X s (M) —> X S ~ 1 (M) is not 
an operation on X s (M). In general, the tangent space at an arbitrary 77 E G s is 
T V G S = \U : M —> TM of class H s \ U(m) E If s > 1 + | and X E X S (M ), 

then its global (since M is compact) flow M 3 t 1 —> F t E G s exists and is a (W-curve 
in G s (see, e.g., m Theorem 2.4.2]). The candidate of what should have been the 
Lie group exponential map is exp : T e G s = X s (M) 3 A" 1 —> F\ E G s , where F t is the 
flow of A; however, exp does not cover a neighborhood of the identity and it is not 
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C l . Therefore, all classical proofs in the theory of finite dimensional Lie groups based 
on the exponential map, break down for G s . From now on, we shall always assume 


*>! + !■ 


Since right translation is smooth, each X £ X(M) induces a C°° right invariant 
vector field X R £ £(G S ) on G s , defined by X R (rj) := X o 77. With this notation, we 
have the identity [X R ,Y R ^ (e) = [X, Y], for any X, Y £ X(M). This is the analogue 
of saying that £(M) is the “right Lie algebra” of G s . 

Assume, in addition, that M is connected, oriented, and Riemannian; denote by 
(•, •) the Riemannian metric. Let /i 9 be the Riemannian volume form on M, whose 
expression in local coordinates (x 1 ,... ,x n ) is /i 9 = ^/det (g^dx 1 A • • • A dx n , where 
gij := ( d/dx\ d/dxi) for all i, j = 1,.. .n. Let K : T(TM ) —>• TM be the connector of 
the Levi-Civita connection (with Christoffel symbols T* fc ) defined by the Riemannian 
metric on M; in local coordinates, this intrinsic object, which is a vector bundle map 
K : T(TM ) -3 TM covering the canonical vector bundle projection Tm '■ TM —y M, 
has the expression K(x\ u\ v\ w l ) = (x l , w l + . The connector K satisfies 

hence the identity tm°K = and has the property that the vector bundles ttm '■ 

T(TM) -3 TM and ker K © ker Ttm -3 TM are isomorphic. These two properties 
characterize the connector. Conversely, the connector K determines V in the following 
manner: V xY := K o TY o X for any X, Y £ 3C(M). 

The Riemannian structure on M induces the weak L 2 , or hydrodynamic, metric 
((■,■)), 011 G s given by 



for any 77 £ G s , £ T g G s . This means that the association 77 1 —> ((-,-)) from G s 

to the vector bundle of symmetric covariant two-tensors on G s is smooth but that for 
every 77 £ G s , the map T V G S 3 U v ha ((U v , •)) £ T*G S , where T*G S denotes the linear 
continuous functionals on TG S , is only injective and not, in general, surjective. This 
weak metric is not right invariant (because of the Jacobian appearing in the change of 
variables formula in the integral). 

The usual proof for finite dimensional Lie groups showing the existence of a unique 
Levi-Civita connection associated to a Riemannian metric breaks down, because ((•, •)) 
is weak; the proof would only show uniqueness. However, K° : T{TG S ) -3- TG S 
given by K 0 (Zu v ) := K o Zu v , where Zu £ Tu v {TG s ), is a connector for the vector 
bundle Ta> : TG S -3 G s (since tq‘ 0 K° = tq^ 0 ttg b an d the vector bundles ttg s '■ 
T{TG S ) -3 TG S and ker K° © ker Tt G b -3 TG S are isomorphic). Here, Z Urt £ T Uri (TG s ) 
means that Zu : M -3 T(TM) satisfies ttm (^■c/ 7 ,( m )) £ The covariant 

derivative V° on G s is defined by X Q x y := K° o Ty o X , for any X, y £ £(G S ). This is 
the Levi-Civita connection associated to the weak metric ((•,•)) since it is torsion free 
(V; \y~XyX = [X, (V], for all X, y £ X S (M), where [X, J/] is the Jacobi-Lie bracket of 
vector fields on G s ) and ((•, ■))-compatible (Z[((X, (V))] = {{VzX,y)) + ((X, for 
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all X,y,Z G 3C S (M))] see [20]. Theorem 9.1]. Uniqueness of such a connection follows 
from the weak non-degeneracy of the metric ((•,•)). There is an explicit formula for 
right-invariant covariant derivatives on diffeomorphism groups (see, e.g., [271 page 6]). 
For V° this formula is 

(v° x y) M ■.= k (yfa) o if 1 ) o ,i + (v Y ,r”) o n, 

where X,y G £(G S ), X v X o ri~ 1 } Y v := y o r/" 1 G X s (M), and t hg r) t is a C 1 
curve in G s such that r] 0 = rj and 77* = X{rf)\ this formula is identical to [551 

(3.1)]. Note that each term on the right hand side of this formula is only of class 
TP -1 and, nevertheless, their sum is of class H s because of the abstract definition of 
the covariant derivative on G s . A similar phenomenon occurs with the geodesic spray 
TG S — y T(TG S ) of the weak Riemannian metric; its existence and smoothness for 
boundaryless M was proved in m, i-e., the geodesic spray is in X{TG S ). If M has a 
boundary, this statement is false. 

The discussion above shows that one cannot apply the theorems of Section [3] to the 
infinite dimensional group G s directly. For infinite dimensional problems, they serve 
only as a guideline and direct proofs are needed, which is what we do below. However, 
for each important formula, we shall point out the analogue in the finite dimensional 
abstract setting of Section [3] which inspired the result and we show that Theorem 3.4 
still holds for the model presented here. 

From now on we consider the case M = T 3 . We focus on the following type of SDEs 
on G s , 

( dg u (t, 9) = Yj)=i Hj(g u (t, 9))6W^(t) + u(t, g u (t, 9))dt 
\ 5^(0, 9) = 9, 9e T 3 , 

where Hj G 3C(GU), u(t) := u(t) — \ u G C' 1 ([0, T]; T e G s ) is non-random, 

and dg u (t,9) is the fto differential of g u (t,9) relative to the time variable t. In 
particular, take the constant vector fields H 1 = \/2z/(l, 0,0), H 2 = v / 2z/(0,1,0), 
H 3 = v / 2z/(0, 0,1) on T 3 , where v > 0 is a constant. This is understood in the 
trivialization TT 3 = T 3 x M 3 , so H\ , H 2 , H 3 : T 3 —y M 3 are constant maps. Dehne the 
process g"(t,9) by 


oW(i,0) = VivdW u (t) +u(t,g"(t,8))dt 

s:M = e. 


where W u (t) is a M 3 -valued Brownian motion and dW u {t) is the Ito differential of W u (t) 
with respect to the time variable. 

By standard theory on stochastic flows (see, e.g., [51]), if u is regular enough, 
i.e., u G U 1 ([0, T]; T e G s ) for some s' > 0 large enough, then g”(t,-) G G s for every 
t G [0, T\. From now on, for simplicity, we always assume u to be regular enough. As in 
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[45, Section 6], let U* be some linear space which can be a space of functions, densities, 
or differential forms on T 3 . The action of G s on U* is the pull back map and the action 
of “Lie algebra” T e G s on U* is the Lie derivative. Let ao £ U*, aft ) := ceo#^(f) _1 , 
and aft) : = E[d(£)]. We first study the equation satisfied by aft), given some specific 
choice of ao € U*. 

If we take a 0 = A 0 {9) ■ d 9 := Yl^=i Ao,i(9)d9i to be a C 2 one-form on T 3 , we derive 
the following result (see also an equivalent expression in [[23], equations (32)-(34)). 
Formula (14. 4 p below is the analogue of (I3.32p . derived here by hand for the infinite 
dimensional group G s . 

Proposition 4.1. Let g"(t) be given by (14. ip . Define 

3 

a u (t,9) := (a 0 g"(t, -) _1 ) (9) = ((g"(t, •)” 1 )*ao) (9) := A UJ (t,9)-d9 :='^A i (t,9,uj)d9 i , 

2—1 

where (g”(t, -) -1 )* denotes the pull back map by g"ft, •)~ 1 , and 

3 

aft, 9) : = E [a w (t, 6 1 )] := Aft, 9) ■ d 9 := '^2 Afit, 9)d9{. 

2 — 1 


Then A satisfies the following SPDE, 

3 

dAfit, 9,u>) = — ^ V2udjAi(t, 9, u)dWf{t) 

3 = 1 

3 

9)djAi(t, 9, u) 


3 =1 


+ A fit, 9,u>)diUj(t, 9))dt + uAAfit, 9,u)dt, i = 1,2,3, 


where we use the notation u(t) := (ufit), u 2 (t), u 3 (t)) and dj and A stand for the partial 
derivative and the Laplacian with respect to the space variable 9 of A fit, 9) respectively. 
Equation (14.2p can also be expressed as 


dAfit, 9) 

(4.3) 


y/2fivAfit,9) ■ dWfit) 

(uft, 9) x curlA U) (t, 9) — V(w(i, 9) ■ A fit, 9)2j dt + vA Afit)dt 


{the term dA u {t , 9) above denotes the ltd differential of A^ ft, 9) with respect to the time 
variable) and we have 


(4.4) 


d t A{t , 9) = (uft, 9) x curlA(f, 9) — X7{u(t, 9) ■ Aft, 9))) + vAAft, 9), 
A(0,9) = A o (9). 
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Proof. We use the methods in [185, Lemma 4.1] and [18] Proposition 4.2], It is not hard 
to see that, for the C 2 (note that we assume u to be regular) spatial process Aj(l, 9,cu), 
there exist adapted spatial processes hij(t,9,uj) and Zi(t,9,uj), 1 ^ i,j ^ 3, such that, 

3 

(4.5) dAj(t, 9 , u) = 9 , u)dWf(t) + Zj(t, 9 , u)dt, i — 1,2,3. 

l=i 

We compute below the expressions of hij(t,9,uj ) and Zi(t,9,uj). 

Notice that by the definition of d w (£, 0), 9))* a u (t, 9) = a(O,0) is a constant 

with respect to the time variable. Since 


where the process Vij(t,9,u) := djgf(t,9,u>), we get for each 1 < j < 3, 

3 

(4.6) d (Yl M*, gZ(t, 9),u)V ij (t, 9, w)) =0, 

2=1 

By (14. 5 p and the generalized Ito formula for spatial processes, EH Theorem 3.3.1.], we 
get 


(CM))*«.M) = E 




dAi(t,g"(t,9),uj) = (hij(t,gZ(t,6),u) + V2udjAi(t,g^(t,e),u)j dWfit ) 

l=i 


(4.7) + ^(t,^(l,d),cn) + z/Aii(t,^(f,d),o;) 


X] ( u i ( t i9 v St,0))d j A i (t,g'f(t,9),u) + V^vdjhij (t,g"{t,6),uj\ ) dt. 


l=i 


Using the SDE (14. ip and the theory of the stochastic flows in [511 Theorem 3.3.3.], we 
have, 


3 

(4.8) dV i: j(t, 9 ,uj) — Y d k Uj(t , g%(t, 9))V kj (t, 9, u)dt. 

k =1 

In particular, the martingale part of the above equality vanishes due to the fact that 
the diffusion coefficients in (14.IK are constant. According to (14.71) and (14.8ft . for each 






1 < j < 3, the time It6 differential is 

\i= 1 

3 

= X (hik(t,gZ(t,6),u) + VtodkAi(t,gZ(t,e),u)^V ij (t,6,w)dW%(t) 

i,k=l 

3 / 3 

+ £ EH*. 9u(t,0))dkAi(.t,g”(t,9),uj) + V2ud k h ik (t, g^(t,9),uj)^ 

i= 1 \fc=l 

+ Zi(*, 0„(f, 0), w) + uAAi(t, g u u (t, 0), w) j 9, u)dt 

3 

+ XI ^k(t,g^t,9),u)d i u k (t,g^(t : 9))V ij (t : e,u)dt 

i,k= 1 

Hence from (14.61) . we derive for each 1 < j, m < 3, 

3 

(4-9) ^(him(t,gZ(t,6),u) + Vtod m A i (t,gZ(t,6))')V ij (t,d,u) = 0 

2—1 

and 

X (*<(*> 0«(*» 0 )) + 

i= 1 V 

3 

(4.10) +X («fc(^(M))dfcA(^M),^ + V2vd k h ik (t,g"(t),u) 

+A k {t,gu(t,9),u)diu k (t,gZ(t,0))) ) Vij(t,0,w ) = 0. 


fc=i 


Since {V^(t, 9, is a non-degenerate matrix-valued process (see, e.g., [5T]), 

from (14.9p we deduce that, for each 1 < i, j < 3, 

hij(t,g"(t,9),u) = -y/2vdjAi(t,g"(t,9),u}). 

Noticing that, cu-almost surely, 9 K y g”(t,9) is a diffeomorphism for each fixed t, we 
get 

(4.11) hij(t,9,uj ) = —V2vdjAi(t,9, to), W9 e T 3 , 

which is the expression for hij(t,9,cu). 
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Since {Vij(t, 9, is non-degenerate, by (14.10ft . for each 1 < i < 3, 


Zi(t,g^(t,d),uj) = - vAAi(t,g"(t,6),u) - (u k (t, g"(t, 9))d k Ai(t, g"{t, 9), w)) 


k =1 


-V2ud k h ik {t, g"(t, 9),u) + A k (t , g"(t, 0),u)diU k {t, g"(t, 9 )) 


Taking (14.lip into the above equation and due to the property that, cc-almost surely, 
6 i—^ gu(t, 6) is a diffeomorphism for each fixed t, we obtain the expression for Zi(t, 9 , co), 
namely, 


(4.12) 


Zi(t,9,ui ) =h>AAi(t,9,u) 

3 

- ^ 9)d k Ai(t, 9, u) + A fc (t, 0, u)diU k (t , 0) ) , V0 e T 3 . 

fc=i 


Combining (14.51) . (14.lip , and (14.12p proves (14.2D . We can check that (14.2p is equiva¬ 
lent to (14.3D by direct computation. Taking the expectation of the two sides of equation 
(14.3D and recalling that u(t) is non-random, (14.4ft follows. □ 

Remark 4.2. In Proposition 14.11 U* is taken to be a space of differential forms on T 3 . 
Note that the action of G s on U* is the pull back map and the action of “Lie algebra” 
T e G s on U* is the Lie derivative. Then, for Hi = \/2z/(l, 0, 0), H 2 = \/2z/(0,1,0), 
H 3 = \/2Z'(0, 0,1) and a(t, 9) = A(t, 9) ■ d 9 we have 


3 

v Hi h,=o, Y. a(t)HjHj = uAA(t,9)-d9, 

3 =1 

a(t)u(t) = ( u(t , 9) x curl A(t, 9) — V(w(i, 9) ■ A(t, 0))) • d 9, 

which implies that the second equation in (13.35ft still holds. 0 

Remark 4.3. By the same procedure as in the proof of Proposition l4.ll if a 3 is replaced 
by another term, such as a function or a density, we can still prove the corresponding 
evolution equation for a(t, 9) := E[(cko g£(t, •)) ) (0)]. 

For example, if ao = Po '■ T 3 — y R is a C 2 function, then a(t, 9) satishes the following 
transport equation, see [IS] , 

j d t a(t, 9) = —u(t, 9) ■ Va(t, 9) + vAa(t, 9), 

[ a(0,9) = p 0 (x). 

If ck 0 = D 0 (9)d 3 9 is a density (volume form), write a(t, 9) = D(t , 9)d 3 9. Then D(t , 9) 
satishes the following forward Kolmogorov equation (or Fokker-Plank equation), 

f d t D(t, 0) = -V • ( Du)(t , 9) + uAD(t, 9 ), 

| D(O,e) = D o (0), 
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where V- is the divergence operator. Moreover, if «o = Do(9)d 3 9 is a probability 
measure, let g”(t,9) be the process satisfying 

dg"{t, 9) = \/2udW u (t) + u(t , gfc(t, 9))dt 


whose initial distribution g”(0,9) is D 0 (9)d 3 9. Suppose that for every t G [0,T], the 
distribution of g^(t,9) is of the form D(t : 9)d 3 9: then D(t,9) satisfies (14.131) . 0 


By the analysis in ]3j Section 4.2], for the (infinite dimensional group) G s , and using 
the right-invariant version of definition (12.111) . we have 




dt 


u(t , 9). 


Let «o := (bo(x), B o (0) • dS, D 0 (9)d 3 9), where bo is a C 2 function on T 3 , B o (0) • dS is 
an exact two-form on T 3 , i.e., there is some one-form Aq(9) ■ d 9 such that 


(4.14) B o (0) dS = d(A o (0) -dfl) = J2 (carlA 0 (9)).d9 j Ad9 k , 


and Do{9)d 3 9 is a density on T 3 . We let U* denote the vector space of all such triples 
(b 0 (x).B 0 (9)-dS, D 0 (9)d 3 9). 

As in [33 equation (7.16)], let l : T e G s x U* M. be defined by 

l(u,b,B,D) = j G^\u(0)\ 2 - D(0)e(D(0),b(0)) -^\B(0)A d 3 0, 

where u G T e G s = X s (T 3 ) is the Eulerian (spatial) velocity of the fluid, b G C 2 (T 3 ) is 
the entropy function, B(0) dS is an exact differential two-form as in (j4. 14[) representing 
the magnetic held in the fluid, D{9)d 3 9 is a density on T 3 representing the mass density 
of the fluid, and the function e(D, b) is the fluid’s specific internal energy. The pressure 
p(D, b ) and the temperature T(D, b) of the fluid are given in terms of a thermodynamic 
equation of state for the specific internal energy e, namely de = — pd(-L) + Tdb = 
D 2 -rjj=jdD + Tdb. It is assumed that c 2 := > 0, where c is the adiabatic sound 

speed. 

Now we proceed as in Section [3] deducing the relevant formulas by hand. As in 
(13. 2 p we define, 


(4.15) J v ° ■““■‘(sSSC.sr.sir) “ E f l(u u (t),b(t),B(t),D(t))dt, 

Jo 

where the stochastic processes gff, 1 < i < 4, are defined by (14.ip with /ij > 0 and the 
same u(t, 9), and the curve 1 1 — y (u u (t), b(t), B(f), D{t)) in U* is dehned by 


(4.16) 


/, n\ rji D gtT (tfl) 

U u(t,9) Tgn^Rgn^o)-! Jt , 

b(t,e) :=E[((g»(t,-)- l Y b 0 ) (»)], 

B(t,9) ■ dS E [(( 9 »(i, V 1 )* (B„(0) ■ dS)) («)] , 
D(t,0)d 3 0 := E [(( 9 "(t, .)-')• (£>o(9)d s 9)) («)] , 
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where (g^(t, -) -1 )* denotes the pull back map by g^(t, -) _1 : T 3 —>■ T 3 . 

For each non-random v G C 1 ([0,1]; T e G s ) with u(0) = v(T) = 0 which is regular 
enough, the perturbation (13.4ft (for right invariant systems) in the formulation here is 
determined by the following flow e BtV (t, •) G G s , see e.g., [3], [IT], 

f = £^v(t,e etV {t,9)) 

\ e £) „(0 ,0)=9 

We say that (gfc 1 , gff, gff, g^ff) is a critical point of J v °’" 0 ’* if for each v G C 1 ([0, T]; T e G s ) 
with n(0) = v(T) = 0, 


rt M2 3 /^M4 1 

u,j \ye,vi ye,v> ye,v> ye,v) 


de 


= 0, 


£=0 


where gg v (t,u>){6) := e E>v (t, g^(t,9)). 

Then we have the following result. 

Theorem 4.4. The semimartingale (gff, gff, gff, gff) is a critical point of 3 v ° ,a °’ 1 if 
and only if the following equations hold for ( u(t ), b(t), B (t), D(t )), 


( 


(4.17) { 




dtU + u-Vu+% + Bxc / ) irlB 
dtb = —u ■ V6 + ^ 2 A& 
d,B — curl(w x B) + 1 U 3 AB 
d t D = —V • (Du) + /x 4 A D 

V • B = 0, 


Ah Am + (/ii - + 2 /ii(VlogZD, Vn) 


where p is the pressure term. 


Proof. By the analysis in [3] Section 4.2], properties (I3.20p and (I3.22p still hold for the 
diffeomorphism group G s . Since the right action of G s here is the pull back map, it is 
not different to check the validity of (I3.23P for our present formulation. Therefore we 
can repeat the argument to show that (gff 1 , gff 1 gf, 3 , gff) is a critical point of J v(,, “ 0,z if 
and only if the following equation holds, 


/ A -1 n\ d SI ^ SI SI 1 / \ 

(418) = + »«*(*) 


SB 


oB(t) 


si 

SD 


oD(t) 


K 


where b(t), B(t), D(t) are defined by (14.161) . and the operator K : T e G s —> T e G s is 
defined by (I3.10p . and where we also use the property S7° H .H, = 0; hence u(t) = u(t). 
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Note that B o (0) • dS = d(-A o (0) • d0) for some one-form A 0 (9) ■ d9, and 


where 


B (t,9) • dS : = E [((g^(t, .) _1 )*(BoW ■ &dS))(0) 
= E[((g«(t,.)- 1 )* d (i4o(fl)- d fl))((?) 

= E[d((^(t,.r 1 )*(^)- d «))(<?) 
= d(A(t,e)-dB)(e), 


A(t,9)-d9:=E ((g' i3 (t,.)- 1 Y(A 0 (e)-de))(e) , curL4(f) = B(t) 


By Proposition 14.11 equation (14.4ft holds for A(t) with viscosity constant v = and 
hence B(t) = curL4(t) satisfies, 


<9 t B = curl(w x B) + /i 3 AB. 


We also have V • B(t) = V ■ (curL4(f)) = 0. 

Moreover, from Remark 14.31 we know that (b(t), B(i), D(t)) satisfies the last four 
equations in (14.171) . Since Hi = y/2pi(l, 0, 0), H 2 = y / 2/ii(0,1, 0), H 3 = y / 2/ii(0, 0,1), 
by the analysis in [3], Proposition 4.1], for every v G T e G s , 


1 x A 

2 5Z ( Vad vHjHj + V# .(ad. v Hj)) = -piAv] 
3 = 1 


therefore we have AT(J^) = —piA(Du). 

From the computation in [45, Section 7], particularly (7.4), (7.18) and (7.19), we 
get, 


51 51 51 51 

ad*-1-o b -|-o B H-o D 

5u 5b SB 5D 

= —Du ■ Vu — uV • (Du) — B x curlB 


Vp 


for some function p, the pressure. 

Combining all the above equalities into (14.181) . we obtain, 


(4.19) d t (Du ) = —Du ■ Vu — uV ■ (Du) — B x curlB — Vp + piA(Du), 


and, since 


d t D = —V ■ (Du) + /i 4 A D, 


we derive the first equation in (14.171) . 


□ 
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Remark 4.5. In the proof of the Theorem we see that the first equation in (13.35ft 
holds (see (14.18)1 ) and, by Remarks 14.21 and [4.31 the second equation is also true. Hence 
Theorem 13.51 still holds for our infinite dimensional model. 


Remark 4.6. The equation satisfied by D is the analogue of the continuity equation 
for mass conservation. We have 


d_ 

dt 


I 


DdO 


I [-V • (Du) + /i 4 A D]dd = 0. 


Remark 4.7. The reason for choosing processes g with different constants /q is that 
the viscosity constants in equation (14.17j) are different. 

Furthermore, if we take /q = /i 4 = v for some v > 0, we get, 


(4.20) 


d t u + u ■ Vu + A Vp + T ]3 x curlB = vAu + 2z/(V log D, Vu) 
d t b = —u ■ V6 + /i 2 A b 
< d t B = curl(u x B) + /i 3 AB 
d t D = -V • (Du) + vAD 
V ■ B = 0 


where u is the fluid velocity, b is the entropy, B is the magnetic held, p is the pressure, D 
is the mass density, T is the temperature, and u, /x 3 , and /i 2 are constants representing 
viscosity, resistivity sand diffusivity, respectively. 

This is essentially the equation given in [TQ1 Page 81], where v is the coefficient of 
kinematic viscosity, assumed to be uniform (recall that the kinematic viscosity is the 
ratio of the shear viscosity to the density of the fluid). The term 2z/(V log D, Vu) does 
not appear there. <0 


Remark 4.8. In particular, if we take D(t) = 1, b(t) = 1 for every t G [0, T] in (14.20p . 
we obtain the following viscous MHD equation, see e.g. im 


{ d t u + u ■ Vu + Vp + B x curlB = vAu 
d t B = curl(u x B) + /r 3 AB 
V ■ u = 0, V • B = 0. 0 


Remark 4.9. The reason why the equation for the density (see (14.171) and (14.131) ) 
contains the Laplacian term, and therefore is different from the deterministic case, is 
due to the fact the the underlying paths g^(t) are now random, with Brownian diffusion 
coefficient of intensity given by the viscosity. Thus we have, for all regular functions / 
on T 3 , 



[ f(g”(t, 9))d 3 9 
J T 3 


E 



(( U -V/) + i/A/)( 9 “(i.9)) 4 3 9 . 


which gives precisely the Fokker-Planck equation for D. 
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Remark 4.10. One can check that the energy decays along the solutions of the com¬ 
pressible equations (I4.20p . For this, the presence of the term 2zz(Vlog.D, V«) is im¬ 
portant. <0 


Remark 4.11. The expression Vlog-D in the balance of linear momentum equation 
(the first equation in (14.201) ) is reminiscent of Brenner’s model ([Tj, [8]) for compressible 
fluids which uses two velocities: a mass velocity, that appears only in the continuity 
(mass conservation) equation, and a volume velocity (momentum per unit mass of fluid) 
derived from the motion of particles, that appears, together with the mass velocity, 
in both the balance of linear momentum and the energy conservation equations. The 
fundamental constitutive equation in Brenner’s model (mass, linear momentum, and 
energy conservation) is the statement that the difference between the volume and mass 
velocities is proportional to VlogD, thus leading to a system with only one velocity, 
but having additional terms which involve V log D. It should be emphasized that if the 
fluid is homogeneous and incompressible, the two velocities coincide. We refer to [25] 
for a detailed presentation and an in-depth mathematical analysis of Brenner’s system. 

While the system (I4.20j) involves a single velocity, the fluid part in (I4.20ji seems to 
indicate that the velocity u is Brenner’s volume velocity and that some of the terms 
in his model may be due to stochastic motion of the fluid particles. It would be 
interesting to investigate the relationship of Brenner’s model to the diffeomorphism 
group version of the semidirect product Euler-Poincare equation for stochastic particle 
paths in Theorem 13.51 A different approach, where both Ito calculus with respect to 
forward and backward hltrations describing forward time evolution and its reversal are 
used, can be found in [50] • It also exhibits links to Brenner’s model. <0 
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